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The spectroscopic factor arises from short-range effects in the nuclear wave func- 
tion. On the other hand, cross sections for external capture reactions, such as 
- - - ^Be(p, 7)®B at low energies, depend primarily on the asymptotic normalization of the 

^ ■ bound-state wave function - a long-range property. We investigate the relationship 

, between potential models and the full many-body problem to illustrate how micro- 



scopic substructure effects arise naturally in the relevant transition matrix element 
and can be (in part) accounted for by a spectroscopic factor. 



00 ■ PACS numbers: 25.40.Lw, 26.65. +t, 21.60.-n, 27.20.+ii 

(N 

^ \ The ^Be(p,7)®B reaction at solar energies (-Ecm < 20 keV) plays an important role in 

' the 'solar neutrino puzzle' since the neutrino event rate in the existing chlorine and water 
^ ■ Cerenkov detectors is either dominated by or almost entirely due to the high-energy neutri- 
^ ! nos produced in the subsequent (3 decay of ^B. Direct measurements of the ^Be(p,7)^B rate, 
O ! however, are extremely difficult, since the absolute cross sections diminish exponentially at 
O ■ low energies. Thus, theoretical studies of this process become very valuable. In addition 
to microscopic theories, such as the nuclear shell model or cluster models, potential models 
provide a popular framework for such investigations. In the latter approach, an average po- 
tential is used to generate single-particle wave functions and to calculate various observables; 
microscopic substructure effects are (partly) accounted for through the use of spectroscopic 
^ • factors. This strategy, however, has recently been called into question by Csoto who 
argued that spectroscopic factors should not be included in potential-model calculations of 
\ the '^Be(p, 7)^B cross section. 

■ Csoto's argument is based on the observation that the ''Be(p, 7)^B reaction at low energies 

depends only on the asymptotic normalization of the ^B bound-state wave function [Q, 
whereas the spectroscopic factor arises from the short-ranged properties of the wave function. 
However, the unstated assumption underlying this line of reasoning, namely that the short- 
range and long-range parts of the wave function are independent, is incorrect. For example, 
the bound state is normalized by J d^^r \il){r)\'^ = 1, a global condition which illustrates 
that a change in the wave function at any radius affects the normalization integral and thus 
the value of the wave function everywhere, including in the tail region. Moreover, there 
exists a more subtle connection between the asymptotic normalization and the short-range 
behavior of the wave function based on the bound-state Lippmann-Schwinger equation [Q. 
In fact, the asymptotic normalization can be given in terms of an integral over the interior 
of the nucleus 0. 

The arguments presented in ref. [M, as well as the connection between short-range and 
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long-range parts of the wave function that was outhned in ref. 0, have motivated us to 
reconsider the proper treatment of microscopic nuclear structure effects in potential-model 
descriptions of external capture reactions. In this paper we study how potential models 
are related to the full nuclear many-body problem. We set up a perturbative scheme for 
the many-body system which, in the lowest-order approximation, gives the potential model 
and discuss higher-order corrections coming from the microscopic nuclear structure. The 
spectroscopic factor will be shown to arise naturally in this framework. 

The eight-body wave functions for the ground state and the proton-^Be scattering state 
are written as l^/'f (pi, p2, Ps, p4, Ps, Pe, Pi)) and |V'|(Pi, p2, Ps, p4, Ps, p6, Pr)), respectively. 
Here k is the asymptotic value of the momentum of the scattering particle relative to the 
^Be system, and we have introduced the Jacobi coordinates 

P,^%^-r,^. forj = l,...,8, (1) 

with rg = 0. The dependence of the wave functions on the center-of-mass coordinate, 
R = p8, can be factored out in the form of a plane wave and will be suppressed throughout 
this paper. 

The eigenstates of the seven-body ''Be Hamiltonian, |'?/'"(pi, p2, ps, P4, ps, pe)), form a 
complete set of wave functions when both bound and scattering states are included. Using 
this complete set, the eight-body bound-state wave function for "^B can be written as 

oo 

l^f (P1,P2,P3,P4,P5,P6,P7)) = Yl I^7(P1'P2,P3,P4,P5,P6))0b(P7) , (2) 

n=l 

where the summation symbol indicates a sum over the bound states and an integral over 
the continuum states. The expansion coefficients, 

0b(P7) = (^7(Pl>P2,P3,p4,p5,p6)|^f (P1,P2,P3,P4,P5,P6,P7)) , (3) 

are obtained by inverting eq. (||). Note that (p^i^pj) is the spectroscopic amplitude as defined 
in ref. Its norm is simply the spectroscopic factor, S = J dp7\(j)^^{p7)\'^ and we will show 
below that ^^^(py) can be identified as a one-body wave function. 

The ^B Hamiltonian, Hg, can be written as = Hj — + J2i=i ^d'"* ~ ''"sD) where 

Hj is the Hamiltonian for the ^Be nucleus and /iy = ""^s X]I=i '^j/ Si=i denotes the 
reduced mass of the eight-body system. Since lip^) and \iPq) are eigenstates of Hj and Hg, 
respectively, we have: 

= E-{^|J-\^i) + - ^ + X^^dr. - VslMi) ■ (4) 

Inserting the expansion given in eq. (Q) for \'>p^), we obtain a set of exact, Schrodinger-like, 
coupled equations for the coefficients: 

y2 oo 7 

{Ei - i??)0S(P7) = -I^^IM + E E^d^^ - rM?)^^{P7). (5) 

^'^ m=l i=l 



3 



The normalization condition for the ^^(py) follows from eq. (0): 

oo „ 

E / rfV7|0B(P7)P = l. (6) 

n=l 

The spectroscopic factor is thus 5* = 1 — Yl!^=2 I '^^P7|'/'b(P7)P differs from one because 
of the n 7^ 1 terms in the sum. 

From eq. (|]) we can extract information on the long-range behavior of the ®B ground 
state wave function. The matrix element {ipj\ Yli=i ^(l^i~''"8|)|'^7^); for n oi m correspond- 
ing to a ^Be bound state, has a range on the order of the potential convoluted with the 
bound-state wave function. Thus it falls off rapidly with increasing py. When both ipl^ 
and if)^ describe continuum states, the matrix element has a long-ranged, but infinitesi- 
mal, tail. Solving eq. (|^) outside the range of the potential, we find that the (f)^{pj) de- 
couple and fall off exponentially, (psipr) exp[— -\/2/X7(i?7 — Eg) pi], where pj = \pj\. 
Since {E^ — Eg) is smallest for n = 1, the long-range behavior of the wave func- 
tion is dominated by ^^^(p?) and the associated asymptotic normalization is given by 
An = limp^^oo 0b(P7)P7 exp[+A/2/i7(i?7 — Eg) pj]. Coulomb effects have been ignored in 
these considerations. 

In order to study the relationship between potential models and the full many-body 
problem, we consider the above set of coupled differential equations for arbitrary distances. 
However, we introduce a well-defined approximation scheme. The off-diagonal matrix ele- 
ments Vnm = (^7 I Sl=i V{\fi~f8\)\il^T)y n 171, will be treated as a perturbation. In zeroth 
order we take Vnm = for n 7^ m, that is, only the diagonal terms survive in eq. (|^). The 
equations for the coefficients 4>b{P7) decouple and take the form of Hartree equations for the 
eighth particle moving in the mean field of the other nucleons. The Hartree potential, Vnm, is 
different for each (j)'^{p-j) since the two-body potential is convoluted over different seven-body 
configurations. Assuming that the ground state of ^B is built on the ground state of ^Be, 
only 0j3(p7) can occur in the zeroth-order expansion of that is, |V'^)o = \ip]) 4>b{P7) ■ 

Coefficients (j)^{pj), with n > 1, are associated with excited states of the ^Be system and 
therefore contribute (to lowest order) only to higher-energy states of ^B. Since the *B states 
are orthonormal, it follows directly that the 0^(p7) which correspond to bound states are 
normalized to one and those associated with continuum states are normalized to a delta 
function. 

From the considerations above, it becomes apparent that each </>^(p7) is normalized and 
satisfies a single-particle, Schrodinger-like, wave equation with a mean field of the Hartree 
type. Therefore, the 0^(p7) can be identified with potential-model wave functions. We 
recall that (j)\{p-j) is, by definition, the spectroscopic amplitude, and thus find that the 
potential-model wave function associated with the ground state of '^B yields the lowest- 
order approximation to the spectroscopic amplitude. 

Since, in lowest order, (j)^{pj) = for n > 1, it follows from eq. (H) that there is no 
first-order correction to the coefficient (f)^{pj); the next correction to (f)^^{pj) is of second 
order in the off-diagonal coupling. Nevertheless, coefficients (f)^{p-j) with n > 1 contribute 
to the ®B ground state in first order, thus the norm of 0b (pr) changes from 5 = 1 to 5* < 1. 
The deviation of the spectroscopic factor from one is of second order in the off-diagonal 
matrix elements Vnm- 
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The '^Be-p scattering state \ipi) can be expanded similarly to the bound state: 

oo 

IV'f (Pl, P2, P3, P4, P5, P6, P?)) = 1^7 (Pl' P2> Ps, P4, Ps, P6))4>l{p7) ■ (7) 

n=l 

Proceeding as before, we obtain a set of coupled equations for the expansion coefficients 
05(P7): 

Y2 7 

(^f - ^7)05(P7) = -^<^S(P7) + E E ^(1^^ - ^8|)|V'7"^)C(P7) (8) 

m=l j=l 

A perturbative treatment of the coupling matrix elements Vnm is possible here as well: In 
lowest order the equations decouple and only (psiPr) contributes to I'j/'f). The (psiPr) with 
n > 1 contribute to the scattering state in first order and the next correction to (l)\{pi) is 
of second order in the Vnm- There is no spectroscopic factor associated with the scattering 
wave function, since is normalized by its value at asymptotically large distances, rather 
than by an integral. 

Considering the long-range behavior of iV^f), we infer the following from eq. (^: If the 
incident proton has a kinetic energy less than that of the first excited ^Be state (£'|=0.429 
MeV 0]), (psiPr) becomes oscillatory for large pi, while the (psipr) with n > 1 are exponen- 
tially damped. Hence, for low reaction energies, the n = 1 term dominates the long-range 
behavior of the scattering wave function: \iPq) — >• |V'7)0s'(p7)- At higher energies, 0§(p7) 
with n > 1 display oscillatory behavior and thus a coupled-channels approach to the problem 
may become necessary. 

The dominant contribution to the direct capture process ^Be(p,7)*B arises from an electric 
dipole transition between the ''Be-p scattering state and the '^B ground state . The relevant 
matrix element is given by: 

8 7 

(^8^1 J]e,n|^f) = (^8^1 J^Qp.lV^f) (9) 

1=1 

oo „ 

/ ^'P7 0S*(P7)P70S(P7) 
n=l ^ 

oo 6 « 

+ E E^^-^^l^^n / "^'P^ (P7)C(P7) , (10) 



1=1 i=l 

oo 



n=l i=l 

m— 1 



where the "Jacobi charges", {i = 1, . . . , 8), are defined such that Yl^=i ^i'^'i = Y^i=i ^iPi 
and es equals Q, the total charge of the system. Since the final-state wave function, \ipg), 
is anti-symmetric and the dipole operator, Yl^=i ^i'^i^ is symmetric, it is not necessary to 
anti-symmetrize the initial-state wave function, Without loss of generality, we can 

thus take the eighth nucleon to be the incident particle. Note also that the diagonal matrix 
elements in the last line, Yl^i=i^iPi\^7) ^ vanish since Ylt=i^iPi odd parity and we 
can choose the '^Be wave functions, \ip^), to be states of good parity. 

To lowest order in the perturbative expansion for the wave functions only terms with 
n = m = 1 survive. Thus we obtain: 



Bl 



8 

5^e,r#f) ^ er / d'pr <Pl* (pr) Pr<PUP7) , (H) 
i=i ^ 
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where the Jacobi charge 



er = fir{Sh^-—]^ (12) 



rrii 



accounts for the fact that the photon couples to both the initial nucleus and to the incident 
proton 0. The lowest-order term of the perturbative expansion corresponds to the one-body 
approximation as given, for example, in ref. [Q. 

The first-order corrections to the matrix element of the dipole operator originate entirely 
from the n = 1, m 1 and n ^ 1, m = 1 terms in the second sum of eq. Note 
that Yl^=i^iPi\'^T) J m > 1, contains dipole excitations of the ''Be core, that is, the 
transition operator acts on the internal coordinates of the seven-body system. These first- 
order effects are short-ranged and will be suppressed at low energies, for which the reaction 
takes place predominately in the tail of wave function. 

While there are no first-order corrections to the first sum in eq. (|T0|), there are several 
second-order effects: First, there are n ^ 1 contributions, which will be suppressed at low 
energies since (psiPy) 't>siP7) exponentially damped for large pi. Second, there are 
corrections to both ^^^(pr) and 4>s{.p7) which are due to the off-diagonal coupling. Since 
(l)\{p-j) and (j)\{pi) depend on only one coordinate, these modifications can be accommodated 
by adjusting the one-body potential in a manner similar to the construction of an optical- 
model potential. Third, ^^(py) is no longer normalized to one but to the spectroscopic 
factor, as follows from eq. (|^). In addition, there are second-order corrections to the second 
sum in eq. such as n 7^ 1, m 7^ 1 contributions. These, again, correspond to dipole 

excitations of the ^Be core and are thus short-range effects, which are negligible at small 
energies. 

From the above considerations it follows that the spectroscopic factor should be explicitly 
included in potential-model descriptions of external capture reactions. When the off-diagonal 
matrix elements Vnm are treated perturbatively, the spectroscopic factor arises naturally in 
the normalization of the bound-state wave function at second order and is present at all 
higher orders. For low-energy reactions first-order contributions and second-order effects not 
contained in the spectroscopic factor are either small or may be included by an appropriate 
modification of the one-body potential. At higher energies these first- and second-order 
corrections may become important and it will no longer be sufficient to accommodate many- 
body correlations through the use of spectroscopic factors alone. 

The arguments for including the spectroscopic factor can be taken beyond perturbation 



theory. At low energies, eq. ([TT|) will be valid more generally since the neglected terms are 
associated with short-range effects while the reaction occurs at large distances. Thus we 
have to evaluate an integral which contains two functions of one variable, a scattering wave 
function, (p^^pj), and a bound-state wave function, (p^i^pf). Formally, we can construct 
Schrodinger-like equations for these coefficients using Green functions. To illustrate this, 
we rewrite the coupled differential equations for the expansion coefficients ^^(p?) (eqs. 
and (H)) as follows: 



V2 

f7 ,1 
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P7) + (^^|$^nin-r8|)|^^)0i(pr)-(K-i?7')0i(pr) 

i=l 

= -E(^rlE^(l^^-^8l)l^")C(P7) (13) 



m=2 i=l 
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V2 

P7 , „ 



«(P7) + E(^"iE^(i^^-^8i)i^r)C(P7)-(K-^7)C(P7) 



m=2 1=1 
7 



-Wl5^nin-r8|)|^^)0i(p7) 



(14) 



j=i 



Here a = B or S, for the bound and scattering states, respectively. The latter equation is 
valid for n > 1 and can be expressed as: 



(15) 



m=2 



where 



^i:l(P7) 



V2 

-— ^ - (E^ - El] 



j=l 



(16) 



1=1 



Let G be the Green function matrix which solves Y^'^=2 ^^nk iP^)^\^iPT ^ P'l) ~ ^^Pi" P'i)^nm,- 



It then follow that 



lP7j 



/CXD 6 
^V7E^W(P7,P'7)(^7I E^d^^ - ^8l)l^7)0i(p'7) (17) 
fc=2 i=l 



is a solution of eq. ([T^). Inserting this result into the right-hand side of eq. (|T3D yields an 
exact one-body, Schrodinger-like, wave equation for ^^(py). For the bound-state case, that 
is, for a = S, this wave function is normalized to the spectroscopic factor. 

To summarize, we have studied the relationship between potential-model descriptions of 
the capture reaction ^Be(p,7)^B and the full many-body problem. To this end, we set up 
a perturbative scheme for the many-body system and identified the potential-model wave 
function as the zeroth-order approximation to the spectroscopic amplitude. We showed that 
the lowest-order term in the perturbative expansion of the relevant transition matrix element 
corresponds to the standard one-body approximation; the spectroscopic factor was found to 
arise naturally in the normalization of the bound-state wave function at second order. At 
low energies, where first-order effects and additional second-order corrections are small, the 
spectroscopic factor allows one to incorporate the most important many-body correlations 
in potential-model descriptions of external capture reactions and should thus be included 
explicitly in the calculations. At higher energies, where contributions from the interior of 
the nucleus become important, additional microscopic nuclear structure effects arise and 
including the spectroscopic factor will no longer be sufficient. 
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